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CALCULATIONAL ASPECTS OF COLUSIONAL EXCITATION 
OF IONS IN HOT SOLAR PLASMAS 


FOREWORD 

This r^ort is based on remarks made by the author to the Workshop on 
Problems and Prospects in Electron-Positive, Ion Physics held at the Joint 
Institute for Laboratory Astrophysics, University of Colorado, July 1 and 2, 
1974. 
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CALCULATIONAL ASPECTS OF COLLiSlONAL EXCITATION 
OF IONS IN HOT SOLAR PLASMAS 


I. ASTROPHYSICAL APPUCATION 

The first purpose of these remarks is to bring to the attention of the atomic 
scattering community the cross section requirements of physicists studying UV 
and X-ray emission spectra associated with active and flare-produced plasmas 
in the sun's corona. The studies that are of specific concern here are an inte- 
gral part of an experimental satellite program (OSO-5, OSO-7) of the Goddard 
Space Flight Center group (Neupert et ai. ' ), and in the final section of this re- 
port we shall give specific categories of excitations corresponding to prominent 
coronal lines which they observe under solar disturbed conditions. 

In this first section I shall outline briefly the astrophysical background lead- 
ing to the kind of atomic cross section needed and to discuss the shortcomings 
of presently available methods and programs. The second section deals with 
our general approach (the distorted wave approximation) to the calculational 
problem. The third section will consist of tw’o parts: one discusses the effect 
of resonances in a very qualitative way and gives a few words on how their 
effect is included through the multiehannel quantum defect theory. In the second 
part we derive the threshold law for electron-atomic hydrogen excitation in the 
distorted wave approximation. The fact that the approximation is able to re- 
produce that very non-trivial result is as convincing an argument as we can 
adduce of the power of the method. 


The presence of extremely hot material in the sun's corona is indicated 
among other things by the omission of EUV and X-ray lines. We are here con- 
cerned with the coronal counterparts of active centers [T^ ^ (3-4) x lO^ '“K, 

N, - 10® -10**' cm"^] and flare produced plasmas [T,. - (10-40) x lO^ °K, 

Ne 10*** - 10** cm"^?. In all coronal plasmas there is a balance between 
ionization by electron impact and recombination by radiative and dielectronic 
processes. Transitions from ions with closed shells feature prominently (cf. 
Fig. 2, p. 12). 

Most X-ray and UV line emission occurs as a result of electron collisions. 
The analysis of these spectra is usually done in terms of an energy average rate 
coefficient 

C (T) = Jf^(v) vQ (v)dvHln! . 
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where f^Cv) is a Maxwellian electron distribution of velocities v corresponding 
to a ten^)erature T. The ultimate atomic parameter needed is seen to be in the 
cross section Q^. . 

It is the evaluation of (i = initial state, j = final state) then that is the 
chief concern of the atomic calculation. We sketch in a few words what has 
obviously been a long story. Van Regermorter^ has given a formula and co- 
efRt '.ents, based on tiie Bethe approximation, for <H>tically allowed transitions. 

The Coulomb-Born approximation has been applied by Burgess, Hummer, and 
TuUy^ and by Blaha* for a variety of cases. As invaluable as these and many 
other calculations are, we believe they are not really sufficiently accurate for 
the needs of the solar disturbed problem for the following reasons: 

(i) Even though the temperature is high, the mean energy of the electrons 
is generally lower than that of tiie orbital electrons of the highly stripped ions 
of relevance here. Thus a busic criterion for the validity of the Coulomb-Born 
approximation is not fulfilled. 

(ii) For the same reason the consistent (by which I mean an approximation 
better than Born-Oppenheimer) inclusion of the exchange symmetry between the 
scattered and orbital electrons is necessary. 

(iii) Because of the considerable excess of nuclear charge over number of 
electrons, a great number of autoionization states (i.e. resonances), of the 
electron-ion system in the inelastic domain start playing an important role. 

This is one aspect of a non-trivial quantion mechanical phenomenon which shows 
up in electron-neutral scattering at energies closer to impact ionization threshold . 
We shall discuss this more below, but suffice it to say now that the high charge 

of the nucleus enables this phenomenon to be included effectively, analytically, ^ 
and it is our hope to be able to include these important developments in the Goddard 
program. 


n. THE DISTORTED WAVE APPROXIMATION 

In the form that we shall be concerned with it here, the distorted wave 
approximations starts from two equivalent exact expressions for the excitation 
cross section. 


k k 


( 2 . 1 ) 
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are exact wave functions corresponding to initial conditions with outgoing 
ra^al wave scattering (upper sign) parts or final state solutions with ingoing 
radial (lower signs) waves, k. and kf are the initial and final momenta. 4*^ j 
are wave functions corresponding to unperturbed initial (final) states, i.e. they 
are solutions of the Schrodinger equation with Vf (V^ ) absent from the Hamiltonian. 

The distorted wave approximation, as we shall use it, consists of replacing 
the exact solution by an approximate solution which is a one channel function 
corresponding to elastic scattering from the initial or final state at the appropriate 
energy. As work of Shelton et al.^ and McDowell et al.^ has shown, this is an 
excellent approximation for inelastic scattering of electrons from neutral atoms. 

By writing down this form of the distorted wave (assuming here that the target 
state has zero angular momentum). 





Y, (n.) V* ,(2. . .. N T 1) 

0 i target' ' 




( 2 . 2 ) 


one can see that since the distorted wave function is antisymmetrized ((i) thus 
the effects of exchange have dynamical meaning in that they will give 
additional terms in the equation determining u^(r). Secondly, since a sum over 
partial waves is involved, there is no assumption of high energy. (In fact the 
lower the energy the better since fewer partial waves need be included until 
convergence is achieved). In addition to these well-known advantages of the 
distorted wave method there are additional circumstances which benefit the 
highly ionized charged target application. Because the attraction to the nucleus 
so outweights the mutual repulsion between the orbital electrons, a very simple 
approximation of the target wave function, , suffices. (It should be 

realized that an approximation is always required for any many-electron target.) 
For the same reason the polarization of the target by the incoming electron, 
which can be very important for neutral targets,® is negligible here. There is 
an implicit restriction that Z not be so large that relativistic effects become 
dominant; however various calculations^® indicate that relativistic effects are 
not too important for Z < 26 (which is the case of iron, almost the highest Z of 
direct concern in this astrophysical application - cf section IV). 


Ideally the functions u^(r ) would be numerical solutions of the exchange 
approximation equation obtained by substituting (2.2) into a variational principle 
for the scattering phase shift. The general programs however that have so far 
been written® use a scaled Thomas -Fermi potential with an orthogonality con- 
straint to simulate exchange.®- ** The approximation seems adequate although 
somewhat less than what one might optimally desire. However we shall not 
pursue that aspect of the matter any further in this report. 


k 
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ni. OTHER ITEMS 


In this section I shall discuss two items which are relevant to the present 
prcri>lems and approach. The first has to do with the incorporation of resonances, 
whose effect is a major one in these applications; the second derives an analytical 
form of a threshold law in the distorted wave approximation which we have not 
previously seen, and thereby gives further evidence of the cogency of the method. 

It is by now well known that in the scattering of eiec rons from atoms 
resonance can occur. The major class of these resonances occurs just below 
the energies ol the various excited state of the target atom. Although the neutron 
counterparts of these resonances are a classic feature of nuclear physics, it is 
only within the last 12 years that they have had their most naturai and mathemati- 
cal explanation in the work of Feshbach^^ in a way wh<'ch if anything is even more 
suitable to the atomic problem; they are accordingly and appropriately called 
Feshbach resonances. Qualitatively they can be described as states of the electron 
plus target system which are only weakly coupled to energetically allowable open 
channels. Since the easiest way of forming these states is by attaching the In- 
coming electron to an excited state of the target system (which by definition 
means the target has one its electrons in an excited orbital), they are sometimes 
called doubly excited states, (this picture also explains why the resonant energies 
are generally just below* the excited state energies of the target). 

Now exactly the same phenomenon can occur in electron-ion scattering 
except that as the charge of the nucieus increases relative to the number of 
electrons, tiie energies of these autoionization states descend ever farther below 
the parent excited state until they finally descend below even lower excited states 
of the target ion. Thus if one is dealing wdth a highly stripped ion in the region, 
say between the first and second excited states, there will be many classes of 
resonances which will span the whole energy region and they can have an im- 
portant average effect. This is well demonstrated in a calculation of Hershkowitz 
and Seaton^ reproduced in Figure 1 wherein the average of the collision strength. 


=k^(2Sj + 1) (2Lj + 1)Q.^, 


(3.1) 


for e - cm excitation is seen to be augmented by 50% when states attached to 
the P** threshold arc included. (I do not believe that the calculation is intended 
to imply that the cross section drops sharply at that threshold. Presumably 
resonances attached to even more excited states would take over.) The effect 
of such resonances can be expected to be even larger for more highly charged 
targets. 
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(b) 


Figure 1. (a) The enhancement of the average collision 
strength in e-CIII excitation due to resonances vdiose 
parent s^te is lp°. (b) The individual resonances from 
the partial waves corresponding to incoming p wave and 
outgoing s wave of the impacting electron. Dashed curves 
are distorted wave results and solid curves are the full 


close coupling results, 
and Seaton, Hef. 6. 
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Fortunately this effect can be included within the distorted wave formalism. 

The point is that the scattering is dominated by the known (Coulomb) solutions 
of the electron from the ionic core. The many-body effects of the core electrons 
can be incorporated in terms of quantum defects below threshold which can be 
analytically continued into Coulomb phase shifts above threshold. One can par- 
tition the scattering matrix into open and closed channels and, using this connec- 
tion, obtain from open channel distorted wave calculations, the closed channel 
(resonant) effect on other channels vdiich are open. 

The analytic and programmatic aspects of these problem have been the 
subject of much work and development particularly at University College, London 
as exemplified in Refs. 6, 11 and 14. (We are most grateful to Prof. Seaton and 
his colleagues in advance for the opporbinity of availing ourselves of these 
programs. The work in question will be done in collaboration with Dr. A. K. 

Bhatia of our group.) 

I conclude this section with a brief derivation of the threshold law for ex- 
citation of the degenerate (N; ^ = 0,1.. N-1) level from the ground state of hydrogen 
by electron impact. This result, which was first obtained by Gailitis and Damburg^s 
in a brilliant but unfortunately not very transparent analysis, is nontrivial not 
only because the threshold cross section is finite, but because the oscillations 
awry from threshold are (in the complex variable sense) a non-analytic function 
of the excess energy. The fact then that one can get the correct functional form 
from the distorted wave aiq>roximation speaks for itself about the power of the 
method. (This analysis also serves as a nice heuristic derivation of this threshold 
law, which as far as we know has not previously been given.) 

We start with a model of the electron hydrogen interaction which was pre- 
viously introduced in studying the much more difficult impact iorization threshold 
problem®: we replace the electron repulsion r”* by (r^ + r^)"^ 


V 


1 1 


— + 



J_ 1 
r, %, + Tj 


(3.2) 


This corresponds to directions in which the scattered electron (r^ ) is opposite 
to that of the bound electron (vj)* a configuration that dominates tiie long range'** 
interactions, which as we shall show determine the form of the analytical result. 
Within this model the excitation cross section 




(3.3) 



corresponds to some kind of an average of the N't cross section. The (unper- 
turbed) initial state is the s-wave part of a plane wave on the target in its ground 
state: 


«I) = 
1 


sin(kjr) R,(rj) 


(3.4) 


*1 2 

and the exact final state in the distorted wave approximation is 


^ = 
N 




(3.5) 


(fj) refers to r times the N**’ hydrogenic s state of the target (R, being the 
ground state). The interaction (3.2a) can conveniently be approximated by^ 


- — <'■1 < 


VSr ^ 


(3.2b) 
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where 


b, = 3N^ 


(3.6) 


is the dipole moment of the inner electron and nucleus seen by outer electron. 

Substituting (3.5) into usual variational principle (which is equivalent to 
forming J d^r^ R (r^) (H - E) '|<^ = 0), one derives 


- ''(r) + F^^(r) = 0 




With tb<) approxiii'.«tion (3.2b) for V, the exact solution of (3.6) is 


h .(" > = 


f < **0 
^ *’0 


(3.8) 


a pure Coulomb move and a linear combination of dipole solutions: 


lk»,r 

e '* 


F 







- y/r Ji(v^) 

V® 


(3.9a) 





(3.9b) 


r 



J & N are Bessel & Neumann functions and 4>^ is i^iase factor depending on B which 
in turn depends on the matching of function and slope of i?F^ to F^ at r = r^ 

(i.e. the continuity of the logarithmic derivative r) ( r^)). The normalization 
constant t) Is such that it guarantees that Inside solution matches to the outside 
solution, which approaches (the S-wave part of) a plane wave at r-*®; it is given 
by®*'® 




(3.10) 


Substituting F^ and its derivative into 



we find 


(3.12) 


Proceeding now very quickly, we find that the inelastic cross section reduces 
to [we are concerned only with the form thus we set B = 0 in Eq. (3.9b)] 




L. f ° 
Jo 


V j(r)sinkjrdr 


where 




(3.13) 


■f 


R„(rj)V(r. rj)Rj(rj)drjSr 


const. 
[r„ + 2r]2 


( 


■'O 1*0 

is seen to be independent of And now we let kj^ be sufficiently small that 


1 


k^r <<o„ <<- 


in the region of space which gives the dominant contribution to the integrals. 
Under these conditions F^(r) can be replaced by a function in- 

dependent of according to (3.9a), whereas the small argument expansion 

(c>f which we only need the real part) is distinctly not independent of k : 

i ® W N 


J.^(k^r)} a ^os 


(3.15) 


Thus from (3.13) 


0 . ^ 

'«iN 


C. f r-co. - 

Jr, I" \ 2 / J <r„ ♦ 2 rf 


(3.16) 
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The first integral devoted by C is independent of ; an asymptotic series can 
easily be generated for the second integral the leading term of which together 
with C then provides the final result 


+CjCos(a^ lnk„+0)l' 

N 

That is the Gailitis-Damburg result! 


(3.14) 


IV. ASTROPHYSICAL CROSS SECTION NEEDS 

In the context of the astrophysical applications discussed here, i.e. the 
deduction of temperature and density in solsu* disturbed regions from the (UV anc 
X-ray) spectral data, the following are relevant transitions for which excitation 
cross sections are needed. I am indebted to Or. S. Kastner and to Dr. K. Phillips 
of Dr. Neupert's group for this information, and especially to Dr. Phillips for 
astrophysical parts of section I. 

Targets 

Element CON Mg A1 Si S Ar Ca Cr Mn Fe Ni 

Z 6 8 10 12 13 14 16 18 20 24 25 26 28 

Transitions 

One-electron 

Is -* 2s, 2p, 3p, 4p, 5p 

Two-electron 

(Is)^ -* (Is np) ‘P, 3p; (Is n s) 

n = 2, 3, 4, 5 (n > 2 in order to asses cascade contribution) 

Many-ele itron wave length region 

(2p)" - (2p)^-> ni 15A ^ 

(2%^ (2p)" - (2s) (2p)^^> .100 - 200 A 

r = 1, 2, 3, 4, 5 
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